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The regularly sampled, sleep-wake behavior of a mouse can 
be described stochastically as a time series in which a binary 
valued random variable fluctuates between the two possible 
states. It is readily appreciated that such a series is not 
entirely random, but instead exhibits significant local serial 
correlation (LSC) across a range of time intervals (lags). 

This correlation results in part from the organization of the 
alternating episodes of the two states into consolidated "bouts" 
of sleep and wake (1). Serial correlation beyond the mean 
length of the consolidated bouts, however, would be of signifi- 
cant interest and could arise in one of two ways; i) as the 
result of a rhythmic process governing the regular recurrence 
of bouts, or ii) as the result of an arhythmic process in which 
the lengths of successive bouts are serially correlated with 
one another and/or have skewed distributions. Theoretically, 
the two processes should be distinguishable on the basis of 
their resultant serial correlation functions. Practically, 
recognition of the differences in correlation requires time 
series of sufficient length and resolution. 

The majority of work examining temporal distribution in 
mammalian behavioral state data has embraced the hypothesis 
that a short-period "ultradian" rhythm was the source of ob- 
served LSC. In their original characterization of the REM 
and NREM sleep states in human subjects. Dement and Kleitman 
observed a REM/NREM alternation with a period of approximately 
90 minutes (2). The regularity of the REM/NREM alternation 
and observations of apparent local periodicity in other 
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psychological and physiological variables led Kleitman to 
postulate the existence of a unifying ultradian oscillator, a 
"Basic Rest-Activity Cycle" (BRAC) , which transcended the 
sleep/wake boundary and influenced the expression of a variety 
of behaviors in both states (3) . This attractive hypothesis 
led to a number of attempts to demonstrate ultradian rhythms 
in physiology. Most of these attempts, however, have been 
complicated by difficulties in (i) defining variables that can 
be realiably measured across the entire circadian cycle and 
that demonstrate clear ultradian rhythmicity, (ii) obtaining 
data series of sufficient length and resolution to allow defin- 
itive statement about the presence or absence of periodicity, 
and (iii) allowing for serial correlation arising from other 
sources in the analysis of time series data and, especially, 
in the interpretation of spectral and autocorrelation results (4). 

In this paper, we report data obtained from study of the 
temporal sleep-wake distribution of the mouse. Three factors 
make our experimental system particularly well suited to study 
of the temporal organization of sleep-wake state and the possible 
contribution of a BRAC. First, we are currently capable of 
collecting continuous measurements of sleep-wake state sufficient 
in length (up to 280 days) and resolution (sampled every 10 
seconds) to allow rigorous determination of the significance 
and consistency of apparent ultradian and circadian rhythmicity. 
Second, the alternation of sleep and wakefulness has obvious 
intuitive validity as a physiological marker of the rodent BRAC, 
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should one exist. Finally, though the marginal (first-order) 
probabilities of sleep and wakefulness in the mouse demonstrate 
clear circadian periodicity, both states are distributed 
throughout the entire circadian cycle; i.e. the probability 
P(4>) of sleep varies with phase, <J> , but almost never becomes 
one or zero. Thus, if a BRAC continues throughout the circadian 
cycle, its effect on the expression of sleep and wakefulness 
and, in particular, on the correlation structure of those 
variables should be measurable at all phases of the mouse's 
circadian cycle. This fortuitous feature of mouse sleep data 
obviates the use of questionable pairings of parameters from 
wakefulness and sleep on the assumption that they represent 
multiple indicators of the same rhythmic process (e.g. REM 
sleep and daydreaming fantasy (5)). 

Nine male mice ( Mus musculus ) were studied while individually 
housed in light-tight, insulated enclosures illuminated by 
externally controlled, water- j acketed fluorescent lamps. For 
this study, data were collected under two lighting conditions; 
either continuous darkness ("DD"), or 12 hours of light followed 
by 12 hours of darkness ("LD12:12"). Food and water were 
available ad libitum. EEG recording techniques and automated 
state-scoring methods were performed as described elsewhere (6). 
State was determined every 10 seconds as either active wake 
(wheel running), quiet wake, REM sleep or NREM sleep. (For the 
study reported here, active and quiet wake were pooled to 
form "total wake", while REM and NREM sleep were pooled as 
"total sleep.") We selected fourteen sections from the full-length 
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MEAN CORRELATION vs. LAG 



LD 

* 

< 





• CD 
- tP C 
u fd rd 

CC rH cp 

ps 



c 

o 

• rH 

-P 

fd 

I — I 

CD 

p 



o 

C CD 
O C 

•H -H 
-P i— I 

o 
c 

D 



£ 

■P O 

•rH P 
M-i 

x: 

0^ -P 
D G 



O 

JC 

E-« 



P 4-H 
O 



> 
fd 
CD 

rd x: 



CD 

p 

fd 

Cu 

Cu 

fd 

c o 



i — i co 

(d fd 



•H 
P 
0) CD 



c 

fd 

Tj CD 



CO 



-P 0) 
i — I o > 
fd > — i *h 

o a -p 

o o 

p d qj 

G < — l 
r d fd h 
CD O 

-P CD O 
fd i — i 
£ 

-H 
-P 
J) 

•J 



U CD 

rG 

O -P 



fd 

•H 

• ts 

ii fd 

u 

p 

| *rH 

P U 
^ ' 



c g; 

-p 

•H 

£ 



O Cl 



co 
C 
O 
•H 

-P -P 
O U 
CD CD 
CO 4-1 
4-1 

a o) 
a 

x: 

CD o 

g d 

-P co 

co o 
c c 

•H 

fd " 
-P G 
G O 
O *H 
O -P 
•H 
CQ 03 

c 

rH O 
CD U 
G 

fd Cn 
Qj G 

•rH 

CD -P 
< 1 X 

•rH tT 1 

g -h 

£ «“H 



data trains for inclusion in the analysis on the basis of (i) 
minimal interruption and lost data, (ii) relative stationarity 
of dominant circadian period (for DD samples) , and (iii) total 
length between 20 and 30 complete circadian cycles. Eight 
samples were selected from the DD condition, six from LD12:12 
(four mice were recorded in both experimental conditions). The 
data within each section were then "folded" to align the cycles 
along circadian phase (LD12:12 sections were aligned with clock 
time, DD sections were aligned using the onset of wheel-running 
activity as the reference phase) . Data were averaged across 
cycles and the average was plotted as an estimate of probability 
of total wake or total sleep as a function of circadian phase (7). 
Local serial correlation (normalized joint second-order prob- 
ability) was estimated across cycles at each circadian phase (8). 
This treatment was repeated at all lags in the range 0 - 500 
minutes at steps of one minute. 

The graphs in Figure 1 depict mean LSC (averaged across all 
phases of the circadian cycle) for each 20-30 day sample plotted 
as a function of lag. Note that LSC, averaged over the six 
LD12:12 and eight DD sections, is significantly positive (i.e. 
greater than 0.10) for all lags out to approximately 23 minutes. 
Most important, it is also apparent from Figure 1 that mean LSC 
is a monotonically decreasing function of lag; there are no re- 
peated peaks and troughs either in the mean correlation curves 
or in any individual correlation curve that might reflect the 
contribution of a rhythmic process having a stable period in 
this range. 
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For purposes of simplification, only the representative 
correlational results obtained at a lag of 12 minutes were used 
to illustrate the subsequent analysis of phase-dependence of 
LSC (Fig. 2). Figure 2 presents the LSC functions at 12 minute 
lag for wakefulness from a representative DD section. The 
probability function (marginal probability of wakefulness) is 
plotted below the correlation function. Note that the probability 
of sleep should be equal to one minus the probability of wake 
and the correlation functions of the individual halves of such 
a binary process should be identical (9). In actuality, small 
imperfections in the state scoring process (occasional unscored 
or missing data) produce some deviation from the expected results. 
However, these differences are sufficiently insignificant to 
allow presentation of the data for wake as representative of the 
entire process. 

The figure demonstrates that, although mean LSC is signifi- 
cantly positive at the 12 minute lag, the value of the LSC 
function is not consistent across the entire circadian cycle. 
Instead, serial correlation is maximal at phases of the circadian 
cycle corresponding to maximal probability of wakefulness. This 
specific phase dependence of the local correlation structure is 
a consistent feature for all mice observed, though the relationship 
is stronger for sections recorded in DD than for those in LD12:12 
(Table 1), (10). Additionally, we found the relationship of 

correlation and phase to be independent of lag for a wide range 
of lags in the ultradian range (10 seconds to 20 minutes). (As 
mentioned, mean LSC approaches zero beyond lags of 20-25 minutes 
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TABLE 1 



Local serial correlation (LSC) of wakefulness at circadian maxima 
and minima. Each entry represents the average of the 90 LSC 
values in the 90 degree phase windows encompassing the maximum 
and minimum probabilities of wakefulness. A t-test was used to 
evaluate the difference between mean correlation in the two phase 
windows . 



Mouse and 



Condition 


MAX 


P (WAKE) 


MIN 


P (WAKE ) 




ALSC 


t-test 




PHASE 


LSC 


PHASE 


LSC 








LD01 


14° 


.656 


264° 


. 160 




.496 




LD0 2 


2° 


.240 


256° 


.164 




.076 


X/S_=3. 556 


LD0 3 


11° 


. 417 


263° 


. 189 




. 228 


X 


LD04 

LD05 


5° 

5° 


.434 
. 298 


212° 

226° 


. 180 
. 224 




. 254 
.074 


(p 55 . 02 ) 


LD06 


9° 


.338 


245° 


. 121 




.217 














II 

X 
[ — ] 


1 . 345 














x= 


. 224 














V 


. 155 














s_= 


.063 














X 






DD01 


7° 


.626 


241° 


.287 




. 339 




DD02 


6° 


. 242 


266° 


. Ill 




.131 


X/S =5. 302 


DD0 3 


2° 


.453 


252° 


. 060 




. 393 


X 


DDO 4 
DD05 


7° 

2° 


. 398 
. 355 


270° 

252° 


. 253 
. 124 




.145 

.231 


(p > .01) 


DDO 6 


16° 


.418 


270° 


. 071 




. 347 




DDO 7 


2° 


. 299 


270° 


. 162 




. 137 




DDO 8 


13° 


. 176 


260° 


. 078 




. 098 





Ix-1 


. 821 


x= 


.228 


s x = 


.115 


s_= 


. 043 


X 
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and significant correlation no longer exists at any phase of the 
circadian cycle.) 

The consistency of positive correlation across lags and the 
absence of significant periodic peaks and troughs in mean LSC 
argue against a stable ultradian rhythm as the principal source 
of local correlation. Further, at no lag in this 10 second to 
20 minute range was local serial correlation found to be sig- 
nificantly negative for any circadian phase. Thus it is unlikely 
that the observed circadian changes in correlation are the result 
of either imprecise estimation of the ultradian period or an 
ultradian oscillator whose period changes stably with circadian 
phase. Although these findings alone do not preclude correlation 
arising from an ultradian oscillator, they do suggest that such 
a contribution is minor relative to the significant, general 
local serial correlation that was observed. To account for 
all observed features of the correlation structure, the variability 
of the ultradian "oscillator" would be so large as to defy the 
definition of the term. 

A more plausible explanation for the observed correlation 
structure is that the phase of the circadian cycle corresponding 
to peak wakefulness is associated with a simultaneous increase 
in mean bout length for both states and/or an increase in the 
serial correlation of the lengths of successive bouts of a given 
state. That is, the mechanisms controlling the expression of 
sleep and wakefulness might be altered at this phase of the 
circadian cycle such that the lengths of successive bouts of a 
given state are correlated to a greater degree than at other 



10 



circadian phases (11). Such a model could account for both the 
dependence of serial correlation on circadian phase, and the 
occurrence of significant serial correlation at all lags in 
this range. Further, serial correlation arising from such a 
source can result in the appearance of rhythmicity in frequency 
spectra or autocorrelations (12), a phenomenon that may account 
for a number of reports of ultradian rhythmicity in physiology. 

Lastly, it is clear that the absence of evidence for ultradian 
rhythmicity in sleep and wakefulness does not completely preclude 
such periodicity in the REM and NREM states. However preliminary 
results using the same analysis techniques presented here show no 
evidence of ultradian periodicity in REM and NREM states in the 
mouse (13). To remain viable, therefore, it would now seem 
that models proposing ultradian periodicity must either abandon 
the extension to a rest-activity cycle as proposed by Kleitman 
or invoke a phylogenetic origin for the rhythm after the evolution 
of the mouse. Given other similarities of mammalian sleep-wake 
organization, this latter assumption is difficult. 
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6. Electrocortigram is recorded using standard implant 



techniques for the mouse [m. M. Mitler and R. Levine, 
Psychop hy siolog y 1 _, 122, 1970] and preprocessed (amplified 
and filtered) using a Grass Model 6 polygraph. Automated 
sleep state analysis is performed using software developed 
in our laboratory [F. A. Vincent, SQEEK - A c omputer 
algorithm for the real time automatic sleep state scoring 
of mouse electroencephalogram , Ph.D. Dissertation, Stanford 
University, 1977] and running on a PDP 11/34 computer 
(Digital Equipment Corp. , Maynard, MA) . 



7. The data are treated as a binary time series X^ , 

i = 0, 1, ..., where X^ = 1 if the animal is in the 
particular state (e.g. wake) at that point and X.^ = 0 
if he is not. The analysis is then performed by first 
arranging X^ in a two-way layout, modulo L, where L = 
the estimated cycle length t (tau) of the circadian rhythm 
(for animals in LD12:12, t = T, the environmental period 
of 24 hours). Thus X ^ ^ = x (kL+ y) for (phase) 

1 = 0,1, ... (L-l) and k = 0,1,2, ... . As the data are 

assumed to be stationary at any phase or phases in the 
circadian rhythm within the sample (consisting of m cycles 
of period L) , the probability function, P(4>^), can be 
estimated as 



P( V 



m-1 



I 

k=0 



X (k, l) 



/ m . 
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8. The local serial correlation is estimated as 






1/2 



where 



m-1 



P(4) £ ,4> £ + j ) k ^ Q X (k,£) X (k,£+j) 



m 



is an estimate of the joint probability function. 

A 

9. For a binary time series, if P and X., =1 correspond 

1 ( KL / 36 ) 

/\ 

to state 1 and and X^ k ^ =0 correspond to state 0, 

S\ rs 

then we can compute P^C^), P ^£'^£ + j^ an< ^ r^ (£ , j ) as in 

7 and 8 previously. 

It follows that: 
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c) 



r 0 ( £ / j ) = 



p o l *r+ i± ii^VV!iii l 



po <<*>^1 fi- g o ) VW^'Vh+jO] 



1/2 



(i - vv - V^l + Vtl'hwj' 

[(l-P 1 <4, ) ,))(P 1 (* l )(l-J 1 (* ( , +j ))P 1 (* ( , +j )] 1 / 2 



P 1 ( VW - WVW 

VW^'Vhl+jd 



1/2 



= ( l , j ) . 

10. The apparent weakening of phase-dependence of correlation 
observed in LD12:12 sections is not thought to represent a 
true effect of lighting condition. Rather it is believed to 
be the result of reduced variation in the phase alignment of 
these sections resulting in a larger portion of the cycle 
in which the marginal probability of wake is near or equal 
to one (p (sleep) = o) . The correlation behaves erratically 
when P is very near to one and this effect results in a 
reduced average correlation across the 90 degree window 
encompassing maximum probability of wake. For DD sections, 
the increased phase variability of the free running rhythm 
results in fewer points where the marginal probability is very 
close to one and consequently the correlation function is 
better behaved. 
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11. Conclusions about the phase-dependence of bout length hinge 
on the operative definition of the bout itself. With the 
10-second resolution of state measurement used in this 
study, single 10 second epochs often interrupt bouts of 
50-120 minutes in length. Whether this conservative 
definition represents a reasonable interpretation of 
underlying biological processes is open to question. 

12. Serial correlation in time series can give rise to the 
same effects in spectra estimated from a finite sample 

as those obtained from a time series consisting of a periodic 
signal plus noise. This result gave rise to the hunt for 
periodicities in economic time series, a hunt which has 
proved quite futile. See, for example, Granger and Hughes' 
re-examination of the Beveridge time-series of annual wheat 
prices (J. Roy . Sta t . Soc ■ A134 : 413-428 , 1971). See also 
J. T. Enright, J. Theoret. Biol., "The search for rythmicity 
in biological time series", £3, 426-468, 1965. 
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(In preparation). 
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